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An exact calculation of electromagnetic scattering from a perfectly conducting parabolic cylinder 
is employed to compute Casimir forces in several configurations. These include interactions between 
' a parabolic cylinder and a plane, two parabolic cylinders, and a parabolic cylinder and an ordinary 

' cylinder. To elucidate the effect of boundaries, special attention is focused on the "knife-edge" limit 

}_( , in which the parabolic cylinder becomes a half-plane. Geometrical effects are illustrated by con- 

' sidering arbitrary rotations of a parabolic cylinder around its focal axis, and arbitrary translations 

perpendicular to this axis. A quite different geometrical arrangement is explored for the case of an 
ordinary cylinder placed in the interior of a parabolic cylinder. All of these results extend simply 
■ to nonzero temperatures. 

m 

> \ 

3 . I. INTRODUCTION 

The Casimir force, arising from quantum fluctuations of the electromagnetic field in vacuum, is a striking mani- 
festation of quantum field theory at the mesoscopic scale. Casimir's computation of the force between two parallel 
metallic plates ^jl gives the classic demonstration of this phenomenon. Following its experimental confirmation in 
the past decade [21 , however, the Casimir force is now important to the design of microelectromechanical systems 0]. 
, Potential practical applications have motivated the development of large-scale numerical methods to compute Casimir 
' forces for objects of any shape [4 6]. In contrast, the simplest and most commonly used analytic methods for dealing 
0^ , with complex shapes, such as the proximity force approximation (PFA), rely on pairwise summations, limiting their 
' applicability. 

|. Recently we have developed a formalism 0, Q that relates the Casimir interaction among several objects to the 
• scattering of the electromagnetic field from the objects individually. This method decomposes the path integral 
representation of the Casimir energy Q as a log-determinant [To| in terms of a multiple scattering expansion, as was 
done for asymptotic separations in Ref. [ll| . It can also be regarded as a concrete implementation of the perspective 
emphasized by Schwinger that the fluctuations of the electromagnetic field can be traced back to charge and 
current fiuctuations on the objects. (For additional perspectives on the scattering formalism, see also references in 
Q.) This approach allows us to take advantage of the well-developed machinery of scattering theory. In particular, 
] the availability of exact scattering amplitudes for simple objects, such as spheres and cylinders, has made it possible 
■ ■ to compute the Casimir force for two spheres 0, a sphere and a plate multiple cylinders and cases with 
more than two objects [Tsl. flGj. This formalism has also been applied and extended in a number of other situations 

Here we expand on recent work [2lj that showed how to apply these techniques to parabolic cylinders, another 
example where the scattering amplitudes can be computed exactly. The limiting case when the radius of curvature at 
the tip vanishes, so that the parabolic cylinder becomes a semi- infinite plate (a "knife-edge"), provides a particularly 
interesting application of this approach. One can also model the knife-edge as the limit of a wedge of zero opening 
angle [2^ ; the two approaches are rather complementary as the former is most amenable to numerical computation, 
while the latter yields approximate analytic formulae via a multiple reflection expansion (which is useful for other 
sharp geometries, such as the cone [13] )■ Edge geometries have also been considered in Refs. [23l. |2^. 
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The remainder of the manuscript is organized as fohows: The Helmholtz equation in the parabohc cyhnder coor- 
dinate system is reviewed in Sec. [Ill and exact formulae are derived for the scattering of the electromagnetic field 
from a perfectly conducting parabolic cylinder. The techniques of Refs. 0,9: then employed to find the electro- 
magnetic Casimir interaction energy for a variety of situations: We let the parabolic cylinder interact with a plane 
fSec. IIII[) . a second parabolic cylinder (Sec. IIV[) . or an ordinary cylinder (Sec. |Vl). In these calculations we consider 
arbitrary rotations of the parabolic cylinder around its focal axis and arbitrary translations perpendicular to that axis, 
transformations that are particularly useful when considering the "knife-edge" limit. We also position an ordinary 
or parabolic cylinder inside a parabolic cylinder (Sec. IVI| . and incorporate the effects of thermal corrections in all of 
these calculations (Sec. |VlT|. 



II. SCATTERING IN PARABOLIC CYLINDER COORDINATES 



We begin with a review of scattering theory in parabolic cylinder coordinates [2^. Because the system is 
translationally invariant in the z direction and perfectly reflecting, we can decompose the electromagnetic scattering 
problem into two scalar problems, one with Dirichlct and the other with Neumann boundary conditions. Parabolic 
cylinder coordinates are defined by 

x = fiX, y=^(A^-M^) z = z, (1) 

such that for r = \/x^\-y^, 

A ^ sgn{x)y/r + y ^ ^ y/r - y , (2) 

where we have chosen a convention where /i > but A can have either sign, which is explained in more detail below. 
This restricted domain is sufficient to include all points in space. Note that here we must take sgn(O) = 1, not 
sgn(O) = 0. At fixed z, the surfaces of constant /i are confocal parabolas opening upward (toward positive y) and the 
surfaces of constant A are confocal parabolas opening downward (toward negative y), as shown in Fig. [T] The scale 
factors (metric coefficients) are hx = /i^ = vA^ + /?, hz = 1. 




FIG. 1: Coordinate curves in parabolic cylinder coordinates. 



We would like to solve the Helmholtz equation, which in these coordinates takes the form 
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where eventually we will set k = in. For k real, we expect oscillating traveling wave solutions, while for k real, we 
expect exponentially growing and decaying solutions. This equation is amenable to separation of variables: 

'i>{X,fi,z) ^ L{X)M{fi)Z{z) . (4) 

Separation of the z variable is trivial, Z{z) = e'*^^^, leaving 

1 d^L 1 d^M 



L dX^ M dfi^ 



{ki-k')iX'+^i'), (5) 



which gives the separated equations 



dfi^ 

where q is a separation constant. The solutions are 



^-ikl-k')X'L = -qL, 
'^'^ {kl-k')^Hl - qM, (6) 



L{X) = D^{X) and L{X) = D^^^i{iX) , 
M{^i) = D^{ijl) and M(^) = D^^^i[jl) , (7) 



in terms of the parabolic cylinder function D^{u). Here 



A = Ay^-2^ Vfc2 - fc2 = Ay 2 + , 

fl = fi\J-2i^/k^ -kl = ^i^J 2y/kl + k2 , (8) 



and 




- 1 j , which implies -^-1=2 , fc2 ' 

The second solution in each case is obtained by making the replacements q — > — g, A — > iA, and — > — iyU, the 
combination of which leaves the differential equations invariant. For | argu| < 37r/4, Du{u) ~ w'^e"" as |u| — >■ oo. 
So for both A and /x, we have solutions that both grow and decay exponentially as their argument approaches positive 
infinity. Since the Cartesian radial distance is r = (A^ + /x^)/2, these are ordinary exponentials (not Gaussians) when 
expressed in terms of Cartesian coordinates. 

If we send A to — A and /i — > — we return to the same point in space, so a wavefunction defined everywhere in the 
plane must be invariant under this transformation. From the identity 

A.(-A) = {^IfD^iX) + ^^r-+^D^,^,{tX) , (10) 

we can conclude that v must be a non- negative integer to satisfy this requirement. We begin with the case 

L{X)M{ti) = D,{X)D,{ili) , (11) 

for I' = 0, 1, 2, 3, • • • . For these values of the parabolic cylinder functions with real arguments are simple rescalings 
of the solutions to the quantum harmonic oscillator, and thus are given by a Gaussian times a Hcrmite polynomial. 

The combined solution in Eq. (fTT]) then takes the form of a polynomial in fi and A times v^^^^, and thus represents 
a traveling parabolic wave in the +7/ direction. 

These solutions represent "regular" waves, the analogs of solutions in spherical coordinates involving spherical 
Bessel functions and spherical harmonics, je{kr)Yjf, {6, (j)). Since we have chosen to restrict fi to positive values, it 
will represent the analog of the radial coordinate r. We will also require "outgoing" solutions to the same differential 
equations, the analogs of solutions in spherical coordinates involving spherical Hankel functions and spherical har- 
monics, h'f\kr)Y^l^{9 , (j)) . As in the spherical case, in the irregular solution the function of the "angular" variable A 
is the same, but the function of the "radial" variable is an independent solution to the same differential equation, 

L{X)M{ti) = D,CX)D^,^i{il) , (12) 
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again for z/ = 0, 1, 2, 3, • • • . Even though they do not blow up at /z = (as the outgoing spherical wavefunctions do at 
r = 0), these solutions are not permissible for /.i < because they are not invariant under the combined substitution 

/i — > —/I and A —A. The solution in Eq. asymptotically approaches a polynomial in /x and A times g^'sA^^'", 
and thus represents an outgoing radial parabolic wave. 
We define the full regular and outgoing solutions 

V;r(r) = e^'=^^D.(A)i?_,_i(M), (13) 
using which the free Green's function becomes, for /i > psj.^ 

fih 

G(ri,r2,fc)= / ^^l^V-(r<)>r(r->). (14) 

Here r< (r>) is the point with the smaller (larger) of and ^2, and we have made use of the Wronskian of the two 
independent solutions for each 

W[I?,(M),I?_,_i(m)] (15) 
The decomposition of a plane wave in regular parabolic cylinder functions is (25j 



where 



^ik r ^ikzZ 

~ cos ^ ^ V 2 



1 , ky -\- ikx , , 



and fcy = ^Jk"^ — k"^ — k1 = z-^/ + fc^ + A:^. Here the logarithm defines the arctangent of kx/ky in the appropriate 
quadrant. Note that the expansion in Eq. ([TS]) converges only for ky > Q. since it is built out of parabolic waves that 
propagate upward. 

To determine the T-matrix wc consider Dirichlet or Neumann boundary conditions at ^ = /io > 0. In the region 
/i > /.Jo we have the scattering solution 

$(r) = V'^-(r) - ,^ ^-(^^0) w,r(r) (Dirichlet) (18) 

for Dirichlet boundary conditions and 

$(r)=C^(r)-*^+i-^^i(^V'r(r-) (Neumann) (19) 

for Neumann boundary conditions, where prime denotes the derivative of the parabolic cylinder function with respect 
to its argument and = 0, 1, 2, 3, • • • . These wavefunctions correspond to the scattering T-matrix elements Tuk^u'kt ~ 
2TT5{k, - k'JS,,,%^, with 

= (Dirichlet), 

^ JZki!^ (Neumann), (20) 

for the process where an incoming parabolic wave propagating in the +y direction is scattered into an outgoing 
parabolic wave propagating radially. 



^ The factor of {— 1)"^ is incorrectly omitted in Ref. [25II . 
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The solutions we have obtained allowed us to construct the complete free Green's function, the decomposition of a 
plane wave, and the scattering T-matrices, which contain all the information we will need to carry out our calculations. 
However, we note that there also exists a second set of solutions, representing the timc-rcverscd scattering process, 

L{X)M{ji) = Z?_„_i(zA)i?_,_i(/i) , (21) 



for 1/ = 0, 1, 2, 3, • • • , which also are unchanged for A — > —A and n — >■ — /x. For real fc, these solutions go like e ' v*^ 
and so propagate in the ~y direction, in contrast to the solutions in Eq. ([IT]), which go like e'v^^^^ and propagate 
in the +y direction. We also have the the corresponding irregular solutions, 

L{X)M{fi) = D^,_^(i~X)D,{iil) , (22) 



where again = 0, 1, 2, 3, • • • . These solutions go like e~'V'°'~'^'='' and thus correspond to incoming radial parabolic 
waves. 

Since a decomposition of the Green's function like Eq. (|14p typically consists of a sum over all scattering solutions, 
one might wonder why this second set of solutions does not appear there. One can formally extend the sum to include 
all values of v, but for ;/ < the v\ in the denominator becomes a divergent gamma function, so that these terms all 
gives zero contribution. However, we can also write the Green's function solely in terms of these solutions as 

G(Ai,^i,zi; A2,Ai2,Z2,fc) = 

I _J_^^k.i.,-.^) ^ _ i^_._l(^Al)i^_._l(^A2)Z?-.-l(A<)iJ.(^A>) , (23) 

and the decomposition of the plane wave as 

1 -1 A-u-l / A\ (-''-1) 

ih r Ak~z \ 



^^.■r ^ I COt | ) i (* A)Z?_._ i (A) , (24) 

which now converges only for ky < 0, since it consists only of waves propagating downward. Using these solutions, 
we could construct the analogous scattering solutions for Neumann and Dirichlet boundaries, which represent the 
process where an incoming radial parabolic wave is scattered into an outgoing parabolic wave propagating in the —y 
direction. 



III. PARABOLIC CYLINDER OPPOSITE A PLANE 



To calculate the Casimir force for a parabolic cylinder opposite a plane, we will need an appropriate expression for 
the free Green's function in terms of plane waves, and expansions translating between these two bases. For y2 > yi, 
the free Green's function can be written in Cartesian coordinates as 



G(ri,r2,A:) 



^g»fcz(22-2l)_L 



2tt 



,i(k^(x2-xi)+ky{y2-yi)) 



h 



(25) 



i\/ + k"^ + • We equate this expression to the Green's function in Eq. (fT4|) . expand 



where ky = \J k'^ — k"^ — k: 

the plane wave e'*' ''^ in Eq. (|25|) using Eq. (|T6| . make the substitution k^ —k^, and then use the orthogonality of 
the regular parabolic solutions to equate both sides term by term in the sum over i/. The result is an expansion for 
the irregular parabolic solutions in terms of plane waves: 



tan ■ 



kyy/SlT cos I 



(26) 



which is valid for y < and = 0, 1, 2, 3, • • • . We have not found this result in the previous literature, though it is 
hinted at in [28[. The quantity in brackets then defines the conversion matrix between outgoing parabolic cylinder 
functions and plane waves propagating in the —y direction. It allows us to propagate the outgoing waves from the 
parabolic cylinder downward to the plane. We displace the origin of the Cartesian coordinates for the plane from the 
origin of the parabolic cylinder coordinates by a distance d in the y direction, which simply introduces a factor of 
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Due to invariance along the time and z directions, we can make independent computations for each k and fc^;, 
and then integrate over both quantities in the final result for the Casimir energy. In the scattering theory approach, 
the calculation can be formulated in terms of scattering amplitudes by considering fluctuating multipoles 0, [2ll |. 
or equivalently by using a generalized T-operator formalism 18|. In the latter approach, which we adopt here, the 
ingredients we will need are the T-matrix elements, the expansion of the outgoing wave in terms of plane waves, and 
the normalization factors appearing in the Green's functions in Eqs. p4p and (|25p The T-matrix elements for the 



parabolic cylinder are given in Eq. (j20p . and the T-matrix elements for the plane are simply = ±1 for Neumann 
and Dirichlct boundary conditions respectively. Finally, we must include the appropriate normalization factor Q 

(^parabolic / (— l)'^ \ ^ 

— ^— i , where we can read off C'p^"'^'"'"^ — , / —— and Cf '"^^" ~ 4 / from the expressions for the free Green's 

function in Eqs. (HH) and ((25)) . 

We can then write the energy per unit length as 



E _ f°° dn r°° dk^, , [ idk. 



logdet I "-^{-ir''''^'^WkMTZU..kAd)) , (27) 

where the matrix determinant runs over i/' = 0, 1, 2, 3, • • • . Here we have defined the translation matrix 

-1 ( tan ^ 

U.,M - -TT^- 0^^''"" ' (28) 

Vj/!V27r cosf 

and, for convenience in later expressions in which we consider different orientations of the parabolic cylinder, we have 
written the reverse translation matrix as U^'k^idy = (—1)'^ Lly'k^{d), with U^k^{d) — Lli,k^{d). The complete energy 
per unit length is then 

where we have dropped a factor of (— l)'"^^" -'^^ since it does not change the determinant. We sum this result over 
Dirichlet and Neumann boundary conditions to obtain the full electromagnetic result. This can be compared to the 
proximity force approximation. 



2 Vmo 

which is the sum of equal contributions from the Dirichlet and Neumann cases. 
We can make the following simplifications in Eq. (P^ : 

• The integral over kx is zero if + z/' is odd, it is symmetric in ly, v' , and the integrand is even in kx- 



r r dk, ^ I r ^ ^ — TT 

• We can replace / — / by — / qdq, where q~ n'^ + kt 

Jo 27r 2tt 47r Jq 

• We can further simplify the integral 



,,2d 



, 2n 

^2iky d 



(tan I) 
cos^ I 

which appears in Eq. (j29p with n = {v + 1'')/^- Here n is always an integer, since the translation matrix element 



oo S COS"^ f 



/ fc2 

vanishes \i u + u' is odd. Setting m = W 1 H — ^ , we have 

V + ki 



oo 



/„,2.^4(-l)" r~ du'^^^^^e-^-V^^\ (32) 



1 
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which is given in terms of the confluent hypergeometric function of the second kind f/(a, b, x) as 

/„.2d = 2(-l)"e-2V^^'^r(^n+i^ [/(^n+i 0,4dv/^^Tfcf^ = 2^k_2„-i(2dy^M^) , (33) 

where ki{x) is the Bateman k- function [27| . 

We now review some results that were reported previously using this formalism [2l[. To connect back to the 
physical configuration, it is convenient to represent the final Casimir energy in terms of the radius of curvature at the 
tip R ~ IJ.Q and the separation H = d — R/2. At small separations {H/R <^ 1) the proximity force approximation, 
given by 

^2 dx __ rw 



hcL 720 [H + x2 /{2R)f 960^2 

should be valid. The numerical results in Fig. [2] confirm this expectation with a ratio of actual to PFA energy of 
0.9961 at H/R = 0.25 (with R = 1). We note that since the main contribution to PFA is from the proximal parts of 
the two surfaces, the PFA result in Eq. (p4| also applies to a circular cylinder with the same radius R. In the opposite 
limit, R = 0, the parabolic cylinder becomes a half-plane, and we can express the T-matrix in closed form as well: 

= -i^X -u\ cos ^ (Dirichlct, /io = 0), = i^^^ sin ^ (Neumann, /io = 0) , (35) 

V TT 2 V TT 2 

the nonzero elements of which can be summarized compactly as 

r.^ = (36) 

where even v corresponds to Dirichlct boundary conditions and odd v corresponds to Neumann boundary conditions. 
We can thus write the full electromagnetic energy for the half-plane perpendicular to the plane as 



— = — gdglogdet(W-(-l)''k_„_,,_i(2<zi7)) = --^, (37) 



where the Bateman k-function is nonzero only for v + v' even, and we have dropped factors that cancel in the 
determinant. The factor of (— 1)"^ in this expression arises from the T-matrix element for the plane. Numerically, we 
find Cj^ = 0.0067415, which is shown by a dashed line in Fig. [5] 

T his g eometry was studied using the world-line method for a scalar field with Dirichlct boundary conditions in 
Ref. [23|. (The world-line approach requires a large-scale numerical computation, and it is not known how to extend 
this method beyond the case of a scalar with Dirichlct boundary conditions). In our calculation, the Dirichlct 
component of the electromagnetic field makes a contribution = 0.0060485 to our result, in reasonable agreement 



with the value of Cf = 0.00600(2) in Ref. [23|. These results are also in agreement with the calculation in Ref. [24| 

It is straightforward to extend this calculation to the case where the parabolic cylinder (of any radius) is rotated 
by an angle 9c around its focal axis, as shown in Fig. |31 In place of Eq. (pS)) . we have 

1 ftan^)'' 
\fvQ^ cos- 



^.fcjrf,gc)= r-^ - ^Te^^ U,M0c)^U,M-0c)- (38) 



However, now the integral over is not symmetric, and the matrix elements with v + v' odd need not vanish. 

Wc again consider the i? — >■ limit in analyzing this result. From dimensional analysis, the electromagnetic Casimir 
energy at i? = takes the now ^c-depcndent form 

^ - (39) 

where H = d ioi R = Q. Following Ref. [1^, which considers the Casimir energy for a scalar field with Dirichlct 
boundary conditions in this geometry, we plot c{Oc) = cos{9c)C{9c) in Fig. |4l A particularly interesting limit is 
9c — > 7i'/2, as the two plates become parallel. In this case, the leading contribution to the Casimir energy should 
be proportional to the area of the half-plane according to the parallel plate formula, £\\/{hcA) = —c\\/H^ with 
C|| = 7r2/720; plus a subleading correction due to the edge. Multiplying by cos^c" has removed the divergence in 
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FIG. 2: The energy per unit length times , £H^ /(hcL), plotted versus H/R for 6c = and i? = 1 on a log-hnear scale. The 
dashed line gives the R = limit and the solid curve gives the PFA result. 




FIG. 3: The geometry of a tilted parabolic cylinder in front of a plane. 
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FIG. 4: The dependence of the Casimir energy on the tilt angle for a half-plane opposite a plane. The half-plane is a parabolic 
cylinder with R = 0, which is oriented perpendicular to the plane for 6c ~ and parallel to the plane for 6c = tv/2. The left 
panel shows the coefficient c{6c) (see text) as a function of dc, with the exact parallel plate result at dc ~ 7r/2 marked with a 
cross. The inset shows the Dirichlet (circles) and Neumann (squares) contributions to the full electromagnetic result. The right 
panel again shows c{dc), but now in comparison to the proximity force approximation (solid line). Note the large discrepancy 
between the PFA and the exact results as 6c — > 0. 
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C{6c) as 6c — >■ 7r/2. As in Ref. H^, we assume c{6c — >■ 7r/2) = C||/2 + (^^ — ''■/2)c„dgo, although we cannot rule 
out the possibility of additional non-analytic forms, such as logarithmic or other singularities. With this assumption, 
we can estimate the edge correction Cedgc = 0.0009 from the data in Fig. |4l From the inset in Fig. HI we estimate 
the Dirichlet and Neumann contributions to this result to be c^^^ = —0.0025 (in agreement with [23| within our 
error estimates) and cf^^^ = 0.0034 respectively. Because higher partial waves become more important as 9c — > 7r/2, 
reflecting the divergence in C{dc) in this limit, we have used larger values of ly^^^ for 9c near tt/2. In Fig. 2] we also 
show a comparison to the proximity force approximation. The PFA is clearly of no use at = Oj since it simply 
gives zero, while at = 7r/2 the PFA gives the correct energy but incorrectly has zero slope, since it misses the edge 
correction. 

We have found that the edge correction is small in the electromagnetic case, as a result of the near-cancellation 
between the Dirichlet and Neumann contributions. By using Babinet's principle, it is possible to show that this 
suppression of edge effects is a general feature of any thin conductor, arising because the leading term in the multiple 
reflection expansion is identically zero [29j . 



IV. TWO PARABOLIC CYLINDERS 



We next consider the force between two parabolic cylinders opening in opposite directions, as shown in the left 
panel of Fig. [SJ We will consider the generalization to arbitrary orientation below. We need to express the outgoing 
waves from one parabolic cylinder in terms of the regular waves for the other. We let f represent the coordinates of 
the second parabolic cylinder, x — x, y ~ —y — d, and z — z. Using Eq. (|26p for f, we have 



1 



/Stt 



dkx- 



tan ■ 



, cos ^ 

Now we use the expansion of the plane wave, Eq. ()16|) . to obtain 



(40) 



1 



dkx - 



(tan I)' 



iy COS^ f 



(41) 



where d is the interfocal separation. We can then obtain the translation coefficient from the quantity in brackets. We 
thus obtain the Casimir interaction energy 



hcL 



dn 
2^ 



dkz 
2tt 



logdet l^^<^"{d)TS,U."u'{d) 



(42) 
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where T*^ and are the scattering T-matrix elements for the two parabohc cyhnders (which can have different 
radii), the translation matrix elements are given by 



1 



(tan I) 



(43) 



and the determinant runs over v^v' — 0, 1, 2, 3, • • • . Here again we have defined Liy^,{d) = {-1)"+"' Uy^,{d)\ where 
tluy' (rf)^ is the reverse translation matrix. We sum the results for Dirichlet and Neumann boundary conditions to 
obtain the result for electromagnetism. The analogous numerical simplifications apply here as in the case of the plane, 
and we can use Eq. now with d instead of 2d, to express the translation matrix elements in Eq. (|43p in terms of 
the Batcman k-function. The extension to the tilted case is also analogous; now the angle of rotation can be different 
for the two translation matrices, corresponding to different angles of rotation for the two parabolic cylinders. We 
can also introduce a translation in the x-direction dx , in addition to the existing translation d in the y-dircction. For 
rotations 6c and 6c of the two parabolic cylinders and x-translation dx, we have 



Uyv'(,d, 9c, 9c, dx) 
Uvu'{d,9c,9c,dx) 



Uuu'id, —6c 



(tan^)' (tan^)' 



-dx) : 



4>+0c 



4>+0c 



(44) 



where we must have d > 0, but dx can have either sign, representing a translation in either horizontal direction. 

By considering two parabolic cylinders of zero radius, we can study the Casimir interactions of two half-planes, 
as illustrated in Fig. [51 These techniques, together with a multiple reflection expansion, were used in Ref. [soj to 
obtain a variety of results in half-plane geometries. We take 9c = 6c ~ 7r/2, so that we are considering parallel half- 




FIG. 6: Two half-planes tilted by angles 9c and 6c, and displaced by dx and d. 



planes, where positive dx gives the width of the region over which they overlap, while negative dx gives a horizontal 
displacement of the edges away from each other. In this case, Eq. simplifies to 



i^ui^' (d, —, '^I'ix^ — Uuv' (d, —, -^jdx 



1 



dv 



(1 +j„)(^+^'+3)/2 



(45) 



Results are shown in Fig. [71 along with approximations valid in two limiting cases: First, for dx very negative, we 
can ignore the vertical displacement. The configuration is then equivalent to the case of 9c = 9c = 0, which gives 
£/{hcL) = — 0.0020856/(i^ at separation d, which is shown as a dashed line in Fig. [71 Second, for dx large and 
positive, we can take the standard result for parallel plates £\\/{hcL) = —TT'^dx/{72Qd^) plus twice the edge correction 
Sedse/ {ficL) = 0.0009/(i^ found above for a half-plane parallel to a plane, which is shown as a solid line in Fig. [71 
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As these examples illustrate, our description of the two half-planes is redundant: Different parameter choices lead 
to the same physical configuration, a property we have used to check our calculations. The numerical convergence of 
physically equivalent configurations can be quite different, however. For example, in the case oi 9 = 6 = 0, when both 
dx and d increase, the Casimir interaction energy decreases, since the half-planes are becoming further apart. In the 
scattering bases we have chosen, however, in d this effect appears directly through a decaying exponential, while in 
dx it appears through the cancellation of an oscillating integrand. As a result, we need to maintain d > 0, but can 
consider either sign of dx- 



fic L 

d, 




FIG. 7: Electromagnetic Casimir interaction energy per unit length for overlapping planes as a function of horizontal displace- 
ment, in units of the vertical separation d. Solid points are obtained from the exact calculation described in the text. The solid 
line connecting them is a rational function fit to guide the eye. The dashed line gives the energy for the limit where the planes 
are edge-to-edge, while the solid straight line gives the standard parallel plate result for the overlap area, plus edge corrections. 



V. PARABOLIC CYLINDER AND ORDINARY CYLINDER 



We next consider the case of an ordinary cylinder outside a parabolic cylinder, as shown in the right 
panel of Fig. [5l In ordinary cylindrical coordinates, we have regular solutions given in terms of Bessel func- 
tions, e^^'^'e^^^ Jti^Jk"^ — kir), and outgoing solutions given in terms of Hankel functions of the first kind, 
g^k^z ^ite jj(X) f^2 _ k^r), both indexed by angular momentum £. We use the expansion of a plane wave in regu- 
lar ordinary cylindrical wavef unctions, 

= 6*'=^" e'^^e'^'^Jei^k^ - klr) , (46) 



e 

t=-oo 



where cf) is defined as in Eq. ()17p and 6 and r arc the ordinary cylindrical coordinates for r. In these coordinates, the 
free Green's function is given by 

where r< (r>) is the smaller (larger) of ri and r2. 

The T-matrix elements for an ordinary cylinder of radius TZ arc given in terms of Bessel and Hankel functions and 
their modified counterparts by Ttk^i'k' = 2'!TS{kz — k'^)6££'Tf^ , with 



V = \ = ) , V (Dirichlet), 
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FIG. 8: Interior parabolic cylinder geometries: Two parabolic cylinders inside one another (left panel) and an ordinary cylinder 
inside a parabolic cylinder (right panel). 



I e — 



H. 



— I 

2 



(Neumann) , 



(48) 



where prime indicates a derivative with respect to the function's argument. 

For an ordinary cylinder outside a parabolic cylinder with separation d between the center of the ordinary cylinder 
and the focus of the parabolic cylinder, we substitute Eq. (|46| into Eq. (|40)) to obtain 



OO 

E 

£=-oo 



VStt J- 

1 



tan ^ 



ftanlV 



^^e'^^JiWB - klr) 



(49) 



where again the quantity in brackets is the coefficient we need to compute the translation matrix. 
We thus obtain the Casimir interaction energy 



£ 

hcL 



dK 

2^ 



dk^ 
2tt 



logdet 1, 



T 



?=-oo 



U,i{d)T['U,,t{d) 



(50) 



where 



4(-l)^ 1 



■ ftanf)'' 

COS I 



(51) 



and the determinant again runs over i/, i/' = 0, 1, 2, 3, • • • . As before, for convenience we have defined Uui{d) = 
{-IY+^UM\ where Wj.£(d)^ is the reverse translation matrix. We can again simplify this expression for numerical 
computation by combining the k and kz integrals, and by exploiting symmetries in ^ — >■ —i and k^ — >■ —kr^. General- 
izations to include horizontal translation and tilt of the parabolic cylinder also work in the same way as before, and 
taking the limit in which the radius of the parabolic cylinder goes to zero gives the Casimir energy for a cylindrical 
wire opposite a half-plane. 



VI. INTERIOR GEOMETRIES 



Up to now, we have considered "exterior" geometries in which the objects are outside one another. However, simple 
modifications of these techniques enable us to also consider "interior" geometries using the formalism of Ref. |3l[. 
(Interior geometries were also considered in js^ . using large-scale computation.) There are two changes required for 
this case: The T-matrix of the outside object must be inverted, and we require the translation matrix connecting the 
regular solutions for the different objects, rather than the one connecting outgoing solutions for one object to regular 
solutions for the other. 
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We first consider two parabolic cylinders inside one another. We parameterize the displacement between their foci 
in parabolic cylinder coordinates by Aq and /xg, as shown in the left panel of Fig. [S] Following Ref. [sl], we can 
derive the translation matrix for regular solutions appropriate to the inside problem. Let r' = r + Tq and consider 
the equation e''^ '' = e^'^"^° e^''"^ . Using Eq. we have 



1 ^ 1 f ci>Y „ 1 1 / ci>y 



V'r(Ao,Mo) 



1 1 / d)\ ^ 



2 i/=0 



(52) 



Now we let t = tan — , so cos — = , 

2' 2 VTTt^ 



, and consider the case where \t\ < 1 to obtain 



— V 1 — t/n 1 — ly ■ 

. ./ — n n ^ >, — n 



!/0=0 



Next we write 



if n is odd 



yTTi^ = V a„i" , where a„ = < (-i)"/"r(^) . . 

^ y — r(_i).g.! " ^^'^^ 



T1=0 



Substituting this result into Eq. (|53p and equating powers of i results in 



C'^(A',m')-E 



vo=0 



which yields the coefficient we need from the quantity in brackets. 
The Casimir interaction energy is then 



F r°° riK r°° rik / 1 - \ 

^ = ^logdet(w-(r.^)-'v..,(Ao,A*o)r.9v..,(Ao,A*o)) 



(53) 



(54) 



(55) 



(56) 



where 7^ (7^ ) is the T-matrix for the outer (inner) parabohc cyhnder, 



V,yz.'(Ao,/xo) = Ki/'(-Ao,Aio) = — -au'-u-vo^'Jai^n, ^J^o) , 



1^0=0 



(57) 



and we have dropped normalization factors that cancel in the determinant. 

For an ordinary cylinder inside a parabolic cylinder, as shown in the right panel of Fig. |8l we again let r' = r + rg 
and consider the equation e^^ '^ = e«'= ''og»fc ''^ i^^i ^^ow we expand the left-hand side in parabolic cylinder coordinates 
and the right-hand side in ordinary cylindrical coordinates, to obtain 

1 oo ftanl) oo °° 

_l^V^-«(A',M') = e*'^^'° E e'^''^e''°'°^^o(v/^^~fcl^o)e*'^^^ ^ e^'V^V,(yfc^^r). (58) 



As before, setting t ~ tan — , so that \/l + = 



cos ■ 



it 



and 



-log-j r-, yields 

I 1 — it 



ik 



"'°e'^°'^°Ji,{Vk^ - klra)e'^'''-e'^^Jt{^k^ - klr) . (59) 



I/'— ' ' £— — 00 tQ — — oc 

We take v derivatives with respect to t and then set i = to obtain 



^L°''(A',A*')- E 



E 



1 



dt" \ 71+72 \i-it 



1 + it 



t=o 



(60) 
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where again the quantity in brackets will give the coefficient we need. Using the generalized binomial expansion, we 
obtain 



1 



l + it 



it 



= lA 



min(2|f|.i/) 

E 

n= u mod 2 



(2KI)! 



n\(2\i\-n) 



(^)!r(i-|^| + ^) ' 



(61) 



where the ± is + for £ > and — for ^ < 0, and the sum starts at n = for even and n = 1 for odd, and then in 
both cases goes in steps of 2. For a configuration where the displacement from the focus of the parabolic cylinder to 
the center of the ordinary cylinder is parameterized in ordinary cylindrical coordinates by distance ro and angle 9q 
from the x-axis, the Casimir interaction energy is then 



E 

hcL 



27r 7-00 27r 



logdet W-(r,^) ' J2 ^Mro,Oom°V 



(62) 



where is the T-matrix for the (outer) parabolic cylinder, Tf is the T-matrix for the (inner) ordinary cylinder, 
and 



(63) 



As an example, we consider a thin wire near the focal axis of a parabolic cylinder with parabolic radius To 
leading order in the needle radius TZ, we only need the ^ = Dirichlet T-matrix element, which goes like 1/ log??, for 
TZ small. Keeping only the leading term in 1/ log(7?.//iQ) and also expanding in the displacement from the focus tq, 
we obtain 



£ 
hcL 



3 



15 



32M;5log5 



16Aio log 



T^rosmWo 



Mo 



256fif, log 



■^r^id - 5cos26lo) 



(64) 



Mo 



Here energy is calculated in comparison to the configuration where the ordinary cylinder is placed at x = 0, y = 00. 
In deriving this result, we have assumed that 7?, <C /Xq is small enough that we can drop terms proportional to 
log(^Q-\/K2 + k^) in comparison to terms proportional to log{TZ/fj,l), since for y^K^ _j_ ^2 I/hq the integrand in 
Eq. (|62|) is exponentially suppressed. The first term gives the (negative) Casimir interaction energy per unit length 
when the wire is at the focus, while the second and third terms give the correction as it is moved a small distance away. 
The angular dependence is exactly as we would expect: as the wire moves closer to the vertex axis of the parabolic 
cylinder (Oq = —tt/2), the energy gets more negative; as it moves away from the vertex axis (^o — '^/'^) the energy is 
less negative, and if it moves in a direction perpendicular to the plane of symmetry of the parabolic cylinder (^o = 
or 9q = tt), the energy is unchanged to first order. As a result, unlike the geometric optics calculation considered in 
Ref. [131 1 lisre we do not see any unusual behavior of the Casimir energy at the focus, which is in agreement with the 
results in Ref. [12] ■ In Fig. [S] we illustrate this result for the case where the radius of the ordinary cylinder and its 
displacement are not small. We choose the same radii for the parabolic and ordinary cylinders as in Ref. (32j , and the 
results we obtain are approximately in agreement with what was found there. We cannot make a precise comparison, 
however, because in that work the parabolic cylinder is of finite size and closed at the far end. 



VII. NONZERO TEMPERATURE 



It is straightforward to extend all of these results to temperature T 7^ 0, a subject that has been of significant 

recent interest [1H|- In each calculation, we simply replace the integral |^ by the sum ^^X^i^Lo' '-'^^^ Matsubara 

frequencies k„ = 27rnfcBT/(?ic), where fc^ is Boltzmann's constant and the prime indicates that the n = mode is 

counted with a weight of 1/2 [Sj]. In the classical limit of asymptotically large temperature, only the n = term 

contributes. The numerical calculation is more cumbersome for T 7^ 0, because for T = we could always make the 

f°° dn f°° dk 1 f°° 

substitution / — / — ~fi\/ n"^ + kl) — >■ — / qdqf{q), where q ~ \/ + j since the quantity we integrate 
Jo 27r J_^ l-K ~ Att Jq 

depends only on q. We find it convenient to continue to use the integration variable q, since the quantity we now sum 

and integrate still depends only on this quantity. We therefore carry out the sum and integral via the replacement 



7 rri 00 ' 

n=0 



dk^ 
2tt 



keT 

he 



dk:, 
27r 



dq 



E 

n=l 



tJ 



=/('?) 



(65) 
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FIG. 9: Casimir interaction energy for an ordinary cylinder inside a parabolic cylinder. The left panel shows the geometry, 
with the focus of the parabolic cylinder placed at the origin. We choose the same radii as Ref. [32l, fJ-o ~ V8 and 7?. = 1, so 
that the vertex line of the parabolic cylinder lies at a; = 0, y = —4. The right panel shows the Casimir interaction energy per 
unit length £/{hcL) as a function of x and y, the displacement of the center of the ordinary cylinder from the focus of the 
parabolic cylinder, where the center of the ordinary cylinder lies within the shaded region of the left panel. 



where [xj denotes the greatest integer less than or equal to x. 

As an example, we consider thermal corrections for a conducting half-plane (a parabolic cylinder with i? = 0) 
oriented perpendicular to a conducting plane, at separation H . In the classical limit, only the n = mode contributes 
and we obtain the energy £/L ^ —kp>TCT=nr,/H, with Ct=od = 0.0472. The Dirichlct contribution to this result 
is C^^o^ = 0.0394, in agreement with Ref. [23|. In Fig. [10] the energy for this geometry is shown as a function of 
temperature. For typical separation distances at room temperature, the thermal corrections arc small. 

6H^ 




FIG. 10: The energy per unit length times H'^, £H'^/(hcL), plotted versus kBTH/{hc) for 60=^ and i? = 0. The solid line 
gives the T — >■ cxd limit determined from the lowest Matsubara frequency. For reference, at a separation of H = 350 nm, a 
temperature T = 300 K corresponds to ksT H / {he) « 0.046. 



VIII. CONCLUSIONS 



There are only a limited set of coordinate systems in which the vector Helmholtz equation for electromagnetism 
can be solved exactly. Taking advantage of one of these few cases, we have obtained complete scattering amplitudes 
for a perfectly conducting parabolic cylinder and employed these results to compute Casimir forces. In principle, 
Casimir forces can be computed in configurations involving parabolic cylinders and other shapes for which scattering 
amplitudes are known, as long as we can obtain the translation matrices, which convert expressions of electromagnetic 
waves between different coordinate basis, appropriate to the individual shapes. Following this procedure, we have 
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computed Casimir forces between a parabolic cylinder, a plane, an ordinary cylinder, and a second parabolic cylinder. 
The formalism is versatile enough to treat situations in which one object is enclosed in the interior of a parabolic 
cylinder, and is also easily extended to finite temperatures. 

We focus special attention to the limit when the radius of the parabolic cylinder goes to zero, and it evolves into a 
semi-infinite plate — a knife edge. In this limit we can quantify the contribution of edges to the Casimir force. By 
examining tilted plates, we can consider a broad range of cases involving interacting edges which should be useful 
to the design of microelectromechanical devices. Until recently, the state of art computation of Casimir forces relied 
upon the PFA, which is demonstrably unreliable for a knife edge: A thin metal disk perpendicular to a nearby metal 
surface experiences a Casimir force described by an extension of Eq. (|37|) . while as indicated in Fig. [2l the PFA 
approximation to the energy vanishes as the thickness goes to zero. Based on the full result for perpendicular planes, 
however, we can formulate an "edge PFA," which yields the energy by integrating dS/dL from Eq. ((57)) along the 
edge of the disk. Letting r be the disk radius, in this approximation we obtain 

/v 
{H + r-^/r^- x^)"^dx "^"^") -hcC ^ti ^ r / {2H^) , (66) 

which is valid if the thickness of the disk is small compared to its separation from the plane. (For comparison, note 
that the ordinary PFA for a metal sphere of radius r and a plate is proportional to r / H'^ .) 

A disk may be more experimentally tractable than a plane, since its edge does not need to be maintained parallel 
to the plate. One possibility is a metal film, evaporated onto a substrate that either has low permittivity or can be 
etched away beneath the edge of the deposited film. Micromechanical torsion oscillators, which have already been 
used for Casimir experiments [H], seem readily adaptable for testing Eq. (|39l). Because the overall strength of the 
Casimir effect is weaker for a disk than for a sphere, observing Casimir forces in this geometry will require greater 
sensitivities or shorter separation distances than the sphere-plane case. As the separation gets smaller, however, the 
dominant contributions arise from higher-frequency fiuctuations, and deviations from the perfect conductor limit can 
become important. While the effects of finite conductivity could be captured by an extension of our method, the 
calculation becomes significantly more difficult in this case because the matrix of scattering amplitudes is no longer 
diagonal. 

To estimate the range of important frequencies, we consider R <^ H and 6c = 0. In this case, the integrand 
in Eq. (f37| is strongly peaked around q « 0.3/i?. As a result, by including only values of q up to 2/H, we still 
capture 95% of the full result (and by going up to i/ H we include 99%). This truncation corresponds to a minimum 
"fluctuation wavelength" A„-,i„ — ttH . For the perfect conductor approximation to hold, An,i„ must be large compared 
to the metal's plasma wavelength Ap, so that these fluctuations are well described by assuming perfect reflectivity. 
We also need the thickness of the disk to be small enough compared to H that the deviation from the proximity force 
calculation is evident (see Fig. [2]), but large enough compared to the metal's skin depth 6 that the perfect conductor 
approximation is valid. For a typical metal film, Ap w 130 nm and 5 ~ 25 nm at the relevant wavelengths. For a disk 
of radius r = 100 /im, the present experimental frontier of 0.1 pN sensitivity corresponds to a separation distance 
H sa 350 nm, which then falls within the expected range of validity of our calculation according to these criteria. The 
force could also be enhanced by connecting several identical but well-separated disks. In that case, the same force 
could be measured at a larger separation distance, where our calculation is more accurate. In the case of overlapping 
planes, the correction to the traditional PFA energy is of a similar magnitude to the total force for perpendicular 
planes in the above example, and thus should also be measurable at these separations. We have shown that thermal 
corrections are generally small at room temperature for typical separations, and furthermore our methods allow these 
corrections to be computed precisely. 
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